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ABSTRACT 

This paper answers the question: Which set of 

modes furnishes a higher fidelity math model of 
dynamics of a multibody, deformable spacecraft 
hinges-free or hinges-locked vehicle modes7 Two 
sets of general, discretized, linear equations of 
motion of a spacecraft with an arbitrary number of 
deformable appendages, each articulated directly 
to the core body, are obtained using the above two 
families of modes. By a comparison of these equa- 
tions, ten sets of modal identities are con* 
structed which involve modal momenta coefficients 
and frequencies associated with both classes of 
modes. The sums of infinite series that appear in 
the identities are obtained in terms of mass, and 
first and second moments of inertia of the append- 
ages, core body, and vehicle by using certain 
basic identities concerning appendage modes. 
Applying the above identities to a four-body 
spacecraft, the hinges-locked vehicle modes are 
found to yield a higher fidelity model than 
hinges-free modes, because the latter modes have 
nonconverging modal coef f ic ient s— a characteristic 
proved and illustrated in the paper. 


I. INTRODUCTION 

The use of appendage modes for simulating 
dynamics and control of multibody flexible space- 
craft is widespread, in as much as they are emi- 
nently suitable for both small angle (linear) and 
large angle (nonlinear) dynamics. To win this 
benefit, however, a simulation engineer must 
retain a sufficient number of these modes for each 
appendage so that the simulation program has 
acceptable fidelity. When there are a large number 
of appendages in a spacecraft, and/or an appendage 
has a large mass and moment of inertia relative to 
those of the rigid core body of the spacecraft, 
the total number of appendage modes for a high 
accuracy model may become unacceptably great 
(Reference 1), possibly diminishing the utility of 
the appendage modes for simulation. Furthermore, 
control systems for a multibody spacecraft are 
most easily designed by considering one axis of 
one body at a time, because different bodies 
generally serve different purposes and so the 
control systems' intrinsic features are generally 
quite different. Having designed them so, to 
ensure they all perform as desired in the mutual 
presence and in the presence of flexibility, a 
compact mathematical model of the entire space- 
craft's dynamics is desired so that the control 
designs can be refined fast and economically about 
all axes. For this purpose, the linear, small 
angle models of spacecraft flexible dynamics are 
just right, and so the engineer could beneficially 
employ the vehicle modes of the spacecraft* 
Hughes^ conceived of two families of vehicle modes 
for multibody spacecraft: "hinges-free" and 
"hinges-locked" vehicle modes (although he does 
not use this terminology). By definition, hinges- 


free modes are obtained by leaving all hinges 
free, that is, unlocked and unforced, so that the 
associated natural vehicle modes may contain 
motion of the articulated bodies relative to the 
inboard bodies. Conversely, in the hinges-locked 
modes, the relative motion of the articulated 
bodies is, by definition, zero, and some force or 
torque is applied at the hinges to keep the motion 
so. In Reference 3, these vehicle modes are form- 
ulated, and their zero linear and angular momentum 
properties, the orthogonality conditions, and the 
associated modal momenta coefficients are theo 
rized . 

A critical question whose answer is sought in 
this paper is: Between the hinges-free and hinges- 
locked vehicle modes, which one furnishes a higher 
fidelity dynamic model, retaining the same number 
of modes in the simulation? To this end, a multi- 
body spacecraft is considered in this paper that 
consists of a rigid core body, and N flexible 
appendages, each articulated directly to the core 
body. Three sets of discrete motion equations of 
this spacecraft are obtained from a continuum set 
by using appendage modes, hinges-free vehicle 
modes, and hinges-locked vehicle modes. To compare 
the last two families of modes, modal identities 
are devised that express the sum of contribution 
of all infinite number of modes in terms of first 
and second moments of inertia of the articulated 
bodies, the core body, and the vehicle, following 
Hughes . The analysis is amply illustrated, and 
definitive conclusions are summarized at the end 
of the paper. Although for concreteness, the paper 
considers a multibody spacecraft with Uvel-i 
articulated bodies (the terminology of Ho 3 ), it 
will be clear Chat the conclusions drawn apply to 
a wider range of multibody spacecraft. 

II. FORMULATION OF CONTINUUM EQUATIONS OF MOTION 

Fig* 1 portrays an N+l-body spacecraft that 
consists of a three-axis stabilized core rigid 
body Bq, and deformable bodies E^,...,En’ each 
articulated directly to the core body. The motion 
equations will be developed with respect to the 



Figure 1. An N+l-Body Spacecraft With 
N Articulated, Deformable Appendages 
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reference point 0 in Fig. 1 which is neither the 
mass center ©q of the body Bq, nor the mass 
center® of the entire vehicle V, This generality 
in the formulation is warranted because the 
NASTRAN modal data corresponding to such muitibody 
spacecraft are often with respect to an arbitrary 
reference node 0, and the mass centers are gen* 
erally nodeless empty points. The mass of each 
body is denoted ( p*0 r 1 f • • • ,N) » the mass of all 
N articulated bodies together, m e ; and the mass of 
entire spacecraft, m; dearly, m 3 mg ♦ m fi . The 
first moment of mass of Bq relative to 0 is eg, 
and those for the hinged bodies (j a l,...,N), mea- 
sured from the respective hinges Oj, are 

denoted c., Similar to c (p 3 0 , 1 , , , . ,V) , the 

—j p 

vector r- emanates from 0 and r.from 0*. Vote 

— 0 — j , J 

that the subscript p covers all bodies, while j 
covers only the articulated bodies. The vectors 
b. (j*l,...,N) originating from 0 locate the 
hinges 0* of the hinged bodies Ej. The first 

moment of inertia of the entire spacecraft, then, 
i s 

N 

£ = Eo + l K*i + £oj£j) 4 £o * (1) 

j j j j j J»1 


matrix C _ . transforms 

-Oj 


where the 

j 

vector c. to a Bg-fixed vector, 
vectors . are expressed in the 


the Ej-fixed 
and 


Bg-fixed 


the 
frame. 


Next, Jg denotes the inertia matrix of the body Bq 
about the 

inertia matrix of the hinged body Ej in its own 


reference point 0, while Jj is the 
-ix of the hinged body Ej in it 
frame about the hinge 0 j . The inertia matrix of Ej 
expressed at the reference point 0 in the B n -fixea 


expre 

frame is denoted J'f and 

-j 


= ^Oj -j ~jO 


( m .b x b x ♦ b x 
‘ J-J -J -J 


l£ 0j £j 




( 2 ) 


where (•)* means the 3x3 skew-symmet rix matrix 
associated with the vector (•). The inertia 
matrix J of the entire vehicle at the point 0 will 
then be 


J 3 


l J° 4 j ♦ j° 

4 — j -0 — e 

J 


(3) 


their products can be ignored in the analysis. The 
external forces and torques acting on the 
spacecraft are the force f^(t) and the 

torque £g(t) acting on Bq at 0, and the 

force £j(t) and the torque £j(t) on each Ej at the 

hinge 0-. The latter pair, (f.,£.) , includes a 
J i . J . J 

distributed force £ . . , t )act mg in the domain of 

the body Ej, and if lj(£j»t)is the only surface 

force acting on Ej, then 

f.(t) * . f . ( r . , t ) dA 

-J 'J "J ( 5) 


aj (t) * /j dA 


an elemental area of 


and 


where dA is 

f. A f . With the aid of the Dirac delta 
* j 3 4 E . 

functior! and its derivative, the distributed force 
l.( r.,t) also represents a distributed moment. 
“J 

Regarding the control forces and torques, those 
acting on Bg are included in the quantities 
and whereas, if a control force or torque 

is produced in the interior domain of Ej without 
acting against the core body Bg, then that is 
included in the pair (f.jj£j)l however, if, for 


instance, the torque is produced by an electric 
motor which rests on Bq at 'the interface Oj and 
exerts on E*, then this is considered separately 
and denoted £ Q j(t) (j=l,...,N), for it produces a 
reaction torque ~&g - ( t ) which acts on Bq. The 
total force _f(t) and* torque ^(t) that act on the 
vehicle are 


X = io 


♦ ) c n . f 

<r -Oj-j 


i 0 * l fXojlj ♦ bjCg.f 


J-r 


( 6 ) 


where , 
control 


of course, £(t) 
torque £g.(t) at 


does not 
the interface 


inc Lude 


°r 


the 


The elastic spacecraft under consideration is 
relatively simple; it is straightforward to 
develop its linear, continuum motion equations 
following Hughes^*^^. The equations governing 
the discrete variables (j=l,...,N) 
are: J 


Anticipating our later needs 
matrix Jq . between the bodies 
in the Ej-fixed frame equals 


the cross inertia 
Bg and Ej expressed 


mV, 


\V - c X w_ - y C_.c X Ji. + y f-C^.u. dm 
-0 0 4 -Oj-j-j 4 


= f 


J_ . A J . - 
-Oj 3 -j 


(C h ) x c. x 


(4) 


As for the motion of the spacecraft, its mass 
center is assumed to perform some orbital motion, 
not coupled with its attitude mot ion under 
consideration. To develop motion equations, the 
locaL orbital frame is taken to be an inertial 
frame. The kinetic quantities of interest 

are: V^(t), the percurbat i onal velocity of the 
reference point 0 over the uniform orbital motion 
at time t; wg(t), the inertial angular velocity 
of BgJ £.( t), the angular velocity of each 
articulated body Ej relative to Bg at the 

hinge 0 ; ; and u.(r.,t), the deformation of E- at 
J — J ~ J , V 

the location r. e£.. These quantities are taken to 

be linear, f irst ^ order, infinitesimal, so that 




+ J 






l :&$>> 


. r - 1 u . dm - g 
J-JT -J A 


(7) 


*AoV 




T • 


• J.fl. 

— 0 -j-j 


j . r . u . dm = g n . * 

J j -j-j 


Sj (j 


1 N) 


where an overdot indicates jdLf f erentiat Lon with 
respect to time. To write the motion equation 
governing the deformation £j(£j*t) of the flexible 
body Ej, denote the related linear stiffness 
operator by L.; the body Ej is allowed to be 
anisotropic and/or nonhomogeneous , and its mass 
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density is denoted o.(£.)« The continuum motion 
equation governing the -3 delormat ion Uj is then: 


L.u. ♦ o. {C. n V_ - ( C . n b . ♦ r.) C.-v- 
-J“J J l -jO-0 -jO-j -j — jO— 0 


£j£j * Hj} 3 (j*l,...,N) 


( 8 ) 


The continuum motion equations (7) and (8) are 
discretized in the next section. 

III. DISCRETIZATION OF CONTINUITM 
EQUATIONS OF MOTION 

Three families of modes will be employed in 
this section for discretization: (1) appendage 

modes, (2) hinges-free vehicle modes, and 
(3) hinges-locked vehicle modes. The use of 
appendage modes is standard; they are employed 
here in order to evaluate the infinite sums that 
appear in the hinges-free and hi nges-Locked modal 
identities in Section IV in terms of mass, and 
first and second moments of inertia of the 
appendages, core body, and the vehicle. 

DISCRETIZATION BY APPENDACE MODES 


and y* (t) is the modal input to that mode: 
lo 

T , 

T* (t) - / U* (r.)J.(r.,t) dA (12) 

JO j “jo “J *J “J 

Eqs. (II) are a generalization of Eq. (33) of 
Hughes^, for the former include articulation 
motion of the appendages* Much more complex and 
general equations than Eqs. (11) are available in 
the vast literature on deformable multibody 
dynamics; see, for instance, the works of Ho*, and 
Singh et ai® on spacecraft with arbitrary tree 
topology. Eqs. (11) nevertheless may boast of 
simplicity which is eminently useful while 
designing the control systems for articulated 
bodies. More importantly though, Eqs. (11) are 
derived here because in Section IV they will aid 
in developing modal identities. To facilitate this 
task, Eqs. (11) are abbreviated by using the 
definitions 


JT 4 
2a 

l£oi4* 


J T ft 
■^OA 

f^oi^oi' 

•••» ^on^onI 


Following Hughes 4 , define the modal momenta 
coefficients P. and H. concerning the 

-jo ”J° a 

appendage (cantilever) modes U j ^ (— j ) t ^ 1€ 

articulated body E-: 

-jo * I. d ™ 

H. = j\r*U a (r.) dm (9) 

-jo J-j-JO -j 

where dm = elemental mass. The coefficient Pj CT is 

associated with linear momentum and H. with 

—jo . 

angular momentum of the mode o at the hinge 
point Oj. The modal angular momentum coefficient 
relative to the reference point 0 (Fig. 1) is 
defined as 


H. 6 b* C-.P. ♦ C_ .H . 

-jo s — j —0 j— jo -Oj-jo 


( 10 ) 


Then the continuum equations (7) and (8) 
discretize to 


m YrT c*uL ’ l Cn-c*5. + I C.. J P. Q. 8 f 
-0 0 4 -Oj-j-j 3 “°J a ""J 0 J a “ 

c*V + J#. ♦ £ C..J.1 ♦ HH?Q. * g 

0 - -0 j -Oj-Oj-j j a -jo jo 

x • T • * M 

c C. n V rt ♦ r .C. a ul. + J .3 . ♦ £ H . Q . * 8- + fin- 

-j-jO-O “ 0 j— j 0 — 0 -j-j a -jo jo a j * 0 j 

P T C.A ♦ H 0 <5 ♦ H T 71 . + Q. ♦ Q c Q. * y* 

-jo— j 0—0 -ja -0 -jo-j ^ju jo^jo jo 


(j— i,...,N; o = l,...,®) 


( 11 ) 


where the superscript T indicates transpose of the 
quantity; Q.^(t) is the modal coordinate 

and Q C . is the frequency associated with 
J° a 

the o-th appendage mode U. (r.) of the body E*J 

—jo -j j 


-j “ ^-ji — j2 •••!' 

4 lSoi’1 So-251 

sf Mi h° 2 ...|. 

h a ° T Mill . I 

fij “ IfijX Sj2 -•)- 


—A 


diag [H ; ••• Hj,!. i A ” d i«8 (ij ••• ijJ 


Sj * (‘Iji Qj2 •••!' 2 a * Ifii ••• 2|,1* 

£ 4 ill - ii. 

4 £ 4 [i-£ 1 

£j ft dUg \ a ] { n c ft dug ( 0 *...a 51 


4 1 - w - 4 


T T 

* 1 

N 1 

(13) 


Eqs. (11) then take this concise form: 
T T 

m “ -*^o + 2 A 0 A + 9 a -f 

*i So * 4 a«A * -1 

& ♦ ♦ J A o A * hJ* a • i a ♦ 

* H *“/i + H.S. + 2, ♦ 0 2 9 » Yj 

— A— 0 — A— 0 —A— A ^A — c ~*A —A 


(14) 
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It is interesting to compare Eqs. (14) with 
Eqs* (43) of Hughes^. For even more compaction of 
these equations, the FoLLowing matrices are 
introduced: 



(15) 


where 1^ is a 3x3 identity matrix. Eqs* (14) 
thereby reduce to the following three matrix 
equations: one governing six overalL degrees of 
freedom of the spacecraft, £y(t); the second 
governing n xl vector of reLative angular 
velocities of N articulated bodies; and the third 
governing the ®xl vector of modaL coordinates 
of appendage modes of all articulated bodies. 


T 

«vv”v * *va"a + ^a 3 a - v c) 


(16. a) 


before, because in the case of U-^(£j)» no torque 
acts, by definition, at the hinge 0j to enforce 
the zero deformation and zero rotation condition, 
whereas in the case of Oj Q (£j)? the immobiLe 
support of the appendage enforces that condition. 
Because of the mobile support of the hinge Oj, the 

total motion W. (r.) of E ; in an inertiaL frame is 
-jv -J J 

V-j } • *0v - % + SoA ,X 4o« - j— j~jv 


+ C. .U . (r . ) ( j = l , . . . ,N; v = l , 

-Oj-jv -j 


0 


(18) 


where the first two terms in the right side are 
because of the translation and rotation of the 
core body in the v-th mode, and the third 

term r% . is caused by the relative rotation of 
-j x jv ' . . 

Ej at the free hinge Oj. The motion of the core 

body in v-th mode is given simply by 




^Ov” -0 


(19) 


«vJv * — A-A * j£Jl = *A (t) * * 0 A (t) (16 - b) 

*Sr * ♦ 9 a + ” l!<t> (16 * c) 

Modal identities associated with the modal momenta 

matrices*^ and H,are derived in Section IV. 

—A —A 

DISCRETIZATION BY HINCES-FREE VEHICLE MODES 

In this technique, the continuum equations (7) 
and (8) are discretized all at once. For this 
purpose, the following modal expansion is 
postulated for the variables in Eqs. (7) and Eq. 
(8) (Reference 3): 


Thus a hinges-free vehicle mode W (jr)spans entire 
spacecraft such that 


W (r) 


*0v ( V' 

if I = lo 

w. (r.), 

if r * b. 

-jv -j 

“J 


(j = l 


? 

Following Hughes , the 6 + n^ rigid modes 
of a spacecraft with articulated bodies are 
1, - r X , and -r X ). Not surpr i s ingLy , 
the elastic^ modes W^(r) (v=l T .,.,«) are 
orthogonal to these rigid modes; that is: 



“ 5,(0 ♦ 

z 

*oA U) ' 


“o (t) 

= Qu(t) ♦ 

z 

in n (t), 
-HJv v 


^(t) 

= ^(c) ♦ 

l 

♦ . n (t), 
— jv V 

U 


’ l) = l % 

(l 
) - 

r . )n (t ) , J 
-j v L 

v=I 





<j“l ,2 


, N ) (17) 


where and ©. are the temporal coordinates 

-0 -0 -J 

for the rigid modes of the spacecraft; the total 
number of articuLation degrees of freedom is n a , 
so there are n a + 6 rigid modes in all. Furthermore, 
Rq is the translation of the reference point 0, 
and 0q is the rotation of the spacecraft, 
both in rigid modes} similarly, 0. is the rotation 
of the hinged body Ej relative td^ Bq at the hinge 
Oj (j-l,...,N) in a rigid mode. The quantities 
Zqv* ^0v' and -jv v=l, are v-th 

modal coefficients contributing, respectively, to 
overall discrete motions V^, and £.} and 

n (t) is the associated modaL coordinate. The 

y 

eigenfunction U. (r.) is that part of the hinges- 

free vehicle mdde, "^denoted W (r), which defines 

— v — 

the deformation of body E- in v-th mode. Although 

U. (r.) satisfies the condition of zero dispLace- 
—jv —j 

ment and zero rotation at the hinge Oj, that is, 
Uj v (0j) 3 0 and i 7V v (0j) - 0, it is not the 
same as the a-ch appendage mode Uj^(rJ used 


J* W (£) dm = 0, J £ X W (£) dm = 0, 

V v V v 

f r X W . (r.) dm = 0 (2i.a,b,c) 

J ; -j-jv -J 


where J*y means the entire vehicle is the domain of 
integration. Eqs. (21) can be verified by 
substituting the expansion (17) in the continuum 
equations (7) with zero right sides. Indeed, 
Eqs. (21a, b) state that the Linear and angular 
momentum residing in a v-th hinges-free vehicle 
mode are zero, whereas Eq . (21c) expresses a zero 
momentum-L ike property of the articulated body Ej. 


These properties can be stated alternately by 
defining modal momenta coefficients (p. ,h. ) for 

o. r J v -jv 

each articulated body and (£ ,h ) for all 


articulated bodies collectively: 


p. - J U. (r.) dm h, = f r X U. (r.)dm 

J j-jv-j -JV J j-j-jv-j ^ 

h° > l [b x c n .£. ♦ C n .h. ] 
—v 4 l — j— Oj^-jv “Oj—jv* 


l ~0 \ 


where h^ is defined reLative to the reference 
—v 

point 0. These may be compared with the 

definitions (9) and (10). The “zero momentum 
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properties (21) then transform to: 

m Yn " " 1 C .c fe . ♦ p a 0 

A 0v - ■—Ov j “Oj-j-jv *v 

- iQv - -Ov j -0 J-0 J“J\> -v - 

— j -jfr*Ov — 0 j— jQ^Q v -J~j v — jv “ 

(j-l,...NS v*l ,...,*) (23) 

The eigenvalue problem which governs a hinges-free 
vehicle mode is given by 

L . U . (r.) - a.w 2 C. n W. (r.) v*l, ..♦,“) 

~j“jv -j J v-jO^Jv -j 

(24 ) 

where u is the frequency of v-th mode. The 
orthogonality conditions are obtained by perform- 
ing the operation £ J U^(r^) (•) ^v over the 

eigenvalue problem ( 2^ ) and recalling the proper- 
ties (23). Here, dv is an elemental volume. One 
then obtains 


f W T (r)W (r) dm ■ 6 
J y v - uv 


(25. a) 


Y { r U T U. dm ♦ hT 4> . } + P^Xn + *5 

4 1 J ju— jv — j Vi — J v 1 ^-y^Ov “V “0 v 


J J 


pv 


(25. b) 


l f i. SjAv dra * koiAj^j* " ij^jSjoiovJ ■ 

+ i] y c j0 J j0 i 0v ♦ 4u^j0^ojijv } l 

- !■ 4 u *0v + *0^**0* ' *0 M £ **0v + *0 U - *oJ " ‘ii. 

(25. c) 

y j" uT L.U. (r.)dv = w 2 5 (25. d) 

4 J . — jy — j— jv — j v yv 

where 4 is the Kronecker delta, 
uv 

Utilizing the modal expansion (17), the zero 
momentum modal properties (21) and (23), and the 
orthogonality properties (25), the continuum 
equations are discretized to these decoupled 
equations which separately govern the rigid and 
elastic modes of the spacecraft: 

■ So - z.% - l £oj£-0j - - (26 - a) 

1% + i©o * l " 4 (26,b) 

-j-joK) + + & - Soj + *j U-i,...,*) 

(26. c) 

+ Vv ' *0x4- * + l ±jv ( *0j + £ j ) + Y v (t) 


These equations are abbreviated by recalling 
appropriate definitions from (13) and (15) and by 


*VR 


a 


’ —A 


♦a 


♦ T 

-n 




©1 


T 


T 

A 

*01 

A 

^01 


» *0 * 

T 

*02 

’ 4o * 

* T 

*02 

.Tl 


- - 


- 


— N1 


3 


K 


0 ) * diag 

T 4 r 


( 20 ) 


Here q is a rigid mode vector, whereas in 

(15) is a vector of overall motion of the space- 
craft; the vector (^from (28) and £ A in (13) 

differ likewise. Eqs, (26) now condense to this 
desired compact form: 

-W^VR * ^VA-A ” -V 

«VA%R * -A^A “ *A * *0A 

i + S* 2 ^ = ioi + + ±A ( &0A + *A> + I 

(29) 

DISCRETIZATION BY HINGES-LOCKED VEHICLE MODES 

Since these modes are defined by forcing the 
articulation motion fl. (j*l,...,N) to be zero, 
they are obtained by a modal analysis of the first 
two equations in (7) and Eq. (8) from which 
the £. (j-i,...,N) terms are ignored. The torque 

actually required to keep the hinges locked can be 
evaluated from the third equation in (7) but that 
is not relevant here. The equations for the modal 
anaLysis are therefore: 


• x* 

"*0 " £ *0 


S f C„ . u . dm • 0 
L — 0 J— J 


£ X Xo * i “o * I / j— 0 j + ^Oj-j^-j dm = - 

ijHj * “j %oio - (£jofej + * h 3 > ’ £ 

(30) 

Eqs. (30), in fact, govern the motion of a free 
spacecraft with cantilevered appendages, so the 
sought hinges-locked vehicle modes are the same as 
the unconstrained modes a la Hughes^* The develop- 
ment here parallels that in the previous subsec- 
tion on hinges-free modes. Accordingly, introduce 
the following modal expansion: 


Jo • k ■ I 5 ■ 


(-1.2 -) (26 . d ) So "4b • l sS d S»<‘>. Sj ■ l »]. (J1) 


where y (t) is the scalar input to v-th mode con- 
sidering all articulated bodies collectively: 

T (O & 7 / UT (r )i (r ,t) dA (27) 

j i J J J J 


where the superscript c reminds us that these 

modal quantities pertain to hinges-locked modes, 
c c 

The quantities ^ and are the translation and 
rotation of the core body Bq in a-th mode. Like- 
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wise, the eigenfunction Uj fl (£j) is the defor- 
mation of Ej in the a-th mode, analogous 

to Uj v (£j) in the case of v-th hinges-free mode, 

except that now a force is exerted at the 
hinge Oi to ascertain that U. (0.) = 0 and 

jr r * J 

7 U. (0.) = 0. The total motion of j-th 

J 

appendage relative to the Bg-fixed frame in 

c 

a-th vehicle mode is denoted W. (r.) and it equals 

-ja -j M 

(cf. Eq. (18)] 


(r.) 4 


"J 


*0c 


- (b. + C_ .r . ) <|>: 

-J “Oj-J ~0a 


The a-th mode of the core body, 
other hand, will be 


♦ C^.in (r .) 


-Oj-ja -j 


(32) 


00 the 


v/* 

— 0a ■— 0 


( f ) 4 y 

V Ha 


‘ HAoa 


(33) 


Thus, like Eq. (20), the a-th mode W^(£) uill 

be (r~) or (r.) depending on the domain 

— 0a —0 -ja — j 

under consideration. Orthogonality of these modes 
with the six rigid modes _1 and -r X , similar to 
Eqs. (21. a, 21. b), can be proved easily. To 
express these conditions in terms of hinges-locked 
modal momenta coefficients, define [cf. Eq . (22)]: 


l C . & f u' (r.) dm h' 4 / r x u^ (r ) dm 
£ ja J . -jo -J -ja . -j-ja -J 




—a V — Oj-ja' 

J (34) 


Then, the above mentioned orthogonality is 
[cf. Eq. (23)1: 

m 4a ‘ £ X 4a + 4 ' ^ 


x c ♦ J ♦* * h° - 0 


£ XOa + i Ifla 


(35) 


The eigenvalue problem obeyed by the a-th hinges- 

locked mode W C (r.) is 
“J a 

2 

L.U 1 : (r.) = u C C. U C (r.) (j-1 N) (36) 

-j-ja -j a -jO-ja -j 

where u C is the associated hinges-locked fre- 
quency. The orthogonality conditions between 
a-th and 8-th modes a re 


T T 

f W C W C dm * T / C, A W C fl (r.) dm * 6 . 

j v -o -8 4 -ja-jO-j0 — J a8 


T 

H e C CC / C X C C X C \ 

y U dm - ♦ (-ioa— ^08 ' ±0a^ *08 } 

j 3 

T 

■ »0 a- -08 ‘ 4 iS 


T 

c c 


T 

C X c 


T 

C X c 


y f L.U*? dv 3 a» C 4 (37) 

4 J , -ja-j-jfl i a 8 


J J 


These are a bit more general than Eqs. (62) of 
Hughes^ • 


With the aid of the expansion (31), momental 
properties (35), and orthogonality properties 
(37), the original continuum equations (9) are 
discretized to 


■» «0 -£*§,- l Soi£?3 = iCt) 

£ X «0 * 4) + l 


X •• T •• • r c I-c 

-j-jO-O 0 j j O^H) -j-j £ -ja a a J 

(j=l,...,N) 


T (h" 1 ) U . + n C + oj C " C 3 + + 7 C (t) 

V ja / — j a a a -^aO— — a0 a a ^ 


38) 


Unlike the hinges-free set of discrete 
equations (28), the last two equations in (38) 
involve a new coupling term called "inertial modal 

angular momentum coefficient” defined as 

- Ja 

h*: 1 S /. r X C.U^ (r.) dm - 
-ja j-j-j0-ja “J 

cX ^*/vXrt + (j.C._ — c X C._b X ), C - + J r X U^ dm 
-j-jO^Oa J— j0 -j-jO-j^-aO j -j-ja (3g) 

c o c 

which is different from h. and h defined in 

-ja -a 

(34), The disturbance input y Q (t) Co each 
a-th mode equals (cf, Eq. (29)] 

T 

U' 

j J 

To compact Eqs. (40), introduce 

.T 


Y C (t) 4 l J \j c . (r.)f .(r. f t) dA 
T a 4 J . -ja -j -j f 


(40) 


if * (if if -I. if 4 Ilf - if I 


(41) 


Then, recalling pertinent definitions from (13), 
(15), and (28), Eqs. (38) contract to 

— W “ V R + ^VA— A 


V C> 


cI T **c 


^VA * -A-A * -A ' S A (t > + SoA (t) 

♦ i! c 2 n e (t> * * -c <42) 

The vector n C (t) and y (t), and the matrices 

c c -e “ c 

arC defined like their hinges-free 
companions in (28). 


IV. MODAL IDENTITIES FOR MULTIBODY 
ELASTIC SPACECRAFT 

Our principal concern is to compare hinges-free 
and hinges-locked modes for their accuracy in 
representing articulation motion. To accomplish 
this aim, an equation will be obtained from each 
of the above three sets of discrete equations 
which will be solely in terms of the articulation 
motion and stimuli. These three equations will 
then be compared to yield identities. 
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First, consider the discrete set (16) based on 
appendage modes. By matrix manipulations, the 
following equation governing can be con- 

structed readily: 

(i (% ♦ ^/s z r 1 ™ A ' 9 -' l l ^ n A = 

♦ £'* 2 r l il - «va h v + s H <«> 

where s is the Laplace variable; theoretically, 
SM A is an ® x n fl matrix, and^[ B an <*» x « symmetric 
matrix; 3 is an n fl xn a inertia matrix, and g H is 
the total hinge torque vector: 


i — VA + »A'^ 4 L 


-<P°M - 1 

— A— vv 



- 4 - «VA — VA ' % = 4 A + &0A (44) 

In (44), 1 is ® x ® identity matrix. Thus for an 

equation of 3_ (s),the coefficient of the hinge 

torque ( 9 ) is 


3 J/' 2 rHiL€‘ 1 r l 


(45) 


Anticipating our later needs, now we shall prove 
that the l im s+“ of the matrix [•] in (45) equals 


1 -sj! 9' 1 0 

— —A — ■ — A — 3 — 


(46) 


Applying the matrix inversion lemma to , its 
inverse is found to be: 


-1 0 _T 

%. - i. -T a (t a ° t a 


m r lc r° 

-W 1 -A 


(47) 


On the other hand, owing to the identities 
(D,E,F)" of Hughes 4 

„T 


1»° *»° 


m l 
e— 


c 

— e 


A "l 


(48) 


Also, by definition of in (15), and by virtue 

of Eq. (1) and Eq . (3) 


”w 


m„ 1 - 

0- -0 


^0 


4> 


0 A 
♦ M = 


(49) 


which reduces <11 ^to 


T 

qi _1 =■ I + 0 (50) 

A comparison of amazes. Continuing 

with the proof nevertheless, call upon the basic 
identities (D,E,F)' V of Hughes^ to derive the 
following new identities associated with the 
articulation degrees of freedom: 


a; a, ■ 4 


si a°. • a. 0 


T 0 

11 1 H = j 
-A -A ^A 


^A = - “ ^VA-w4) ^VV -VA, 

^<4 - *} * (o 

(New identities derived in this paper will be 
labeled with Roman numerals as they are cited.) 
These identities and Eq. (50), in turn, lead to 
the identity 

^ ^C^A " - (ID 


which proves Eq . (46) 

An equation analogous to Eq . (43) is obtained 
from the hinges-free discrete set (29). For chat, 
recall the second expansion in (17). Then it can 
be shown that 


i 5 a = ^ » A (i» + '±/* 2 ) 1 (io i + 4) & * l) 

- M M u (s) 

-VA -VV -V v 

* ( I ♦ £ ±J (i. ♦ * 2 / * 2 V l i A l & H (s) (51) 

The coefficient of the hinge torque (s) in 

Eq. (51) equals the term (45). They both reduce 
to 1_ for the lim s+0, and for the lim s+® they 
yield, in view of Eq . (46), the identity 

(i * 0 * A * A r l - 0 dm 

which proves a fortiori that, since the inertia 
matrix $ is positive definite and nonnegative 

definite, the modal coefficients 

v=l, ...,«) [Eq. (17)] constitute a rtonconverging 
series • 


The hinges-locked discrete set 
this aquation for Q^: 

ii ' h* 1 (I * « 2 /s 2 fi A ’ 


(42) furnishes 


M M~ 1 u 
— VA -W -V 


T 

- hi 1 (I ♦ “c/s 2 )' 1 (^ L * a ♦ Ic) (52) 

The equality of the coefficient matrices 
of & H (s) in Eq. (51) and Eq . (52) delivers this 

identity in the s-domain: 

[l - h^ 1 (1 ♦ cii^/s^ )~ 1 h^ I 0~^ ]”^ * l ♦ 

0 *][(.!„,♦ “ 2 /s 2 >" l » A (IV) 

For the lim s-*0, the left side of (IV) degenerates 
to 1 as does its right side. On the other hand, 
caking its limit and recognizing the identity 

(III) produce the identity 

[1 “ * 2 < v > 

The identity (IV) can be rearranged such that 
it reveals poles and zeros of the 
dynamics. For thaj, recognize that when 

i - t ]u (v=l,...,*;j 3 -1) the right side 
of (IV), which is also the coefficient 

of ^(s) in (51), is unbounded, so t jw ^ are the 
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poles of Che spacecraft. Consequently, for 
unboundedness to occur, the Left side of (IV), 
expressed in terms of individual hinges-locked 
modes, yields the identity 

det [l - l (l - u c /w 2 )" 1 h cI h cI >9 _1 1 = 0 (VI) 

* — u v a v —a —a — J 

a 


c X ^ 

where h is the a-th row of the matrix 
c I c 

♦ Similarly, when a ~ ± jui (a = 1,...), the 

matrix within [•] on the # left side of (IV), which 
is the coefficient of 3 is unbounded, 

which impLies that ± are the zeros of the 

dynamics. Therefore, to realize unboundedness, the 
right side of (IV) bears forth 


dec [l ♦ J (l - <J 2 /w C )" 1 5 ♦ * T 1 = 0 (VII) 

4 — v p a ' — — ; ;j — u 1 

T . 

where $ is the u-th row of the matrix 

j> A [Eq. (28)]. Knowing the poles and zeros, the 
identity (IV), keeping in mind its lim, has this 

alternate form (cf. (Y) of Hughes^]: 


t T _ , . n e . 2 

U , cl ,cl n -li-l f / 2^ c x 

■ “ -A -A ^ I * , <» > 

a*l 


, (, 2 +„ 2 ) 

W 


(l * ] « f (s ) 


-A -A 


a=l 


« (s 2 *u 2 ) 

li-l W 


(VIII) 


where 3 total number of retained modes. Because 
of the identity (V), however, this form seems to 
be less useful than the form (IV). Following 
Garg^, one can examine how far the identities (IV) 
or (VIII) are satisfied in the s-domain. The 
identities (VI) and (VII) are useful in several 
ways; for instance, known hinges-locked parameters 
can be used to determine hinges-free modal 
parameters, or vice versa, after Hughes and 
Garg^. Incidentally, the identities (VI) and 
(VII) are analogous to the identities (M) 0 and 
(Q) of Hughes^. As in Reference 4, under 
conditions of symmetry, these identities 
reduce to those concerned with individual 
articulation degrees of freedom. Owing to 
symmetry, <ince different sets of modes will 
contribute to different articulation degrees 
of freedom, the set a (a«l , * . • ,«) may form n^ 
subsets a.(j=l,...,n ) and each a. will span the 
range 1,.A,®; the set m (u 3 1,..*,») fragmentates 
likewise. The identities (VI) and (VII) then 
simplify to 



(X) 


V. ILLUSTRATION OF MODAL IDENTITIES 
AND DISCUSSION 


The identities will now be illustrated for a 
four-body deformable spacecraft shown in Fig. 2. 
It has two flexible soLar arrays, E^ and E 2 , each 
having one articulation degree of freedom about 
yj- and y 2 ~axis, respectively, relative to the 
0-’ _ 

3" 


core body 
degrees of freedom 
four articulation 
0,.., 9^. » and 0 


two rotational 
y^-axis. These 


a sensor having 
about xj- and 
angles are denoted 

**2y f °3x* ® tia °3y* an< ^ the spacecraft thus 
has ten rigid modes. Hinges-free and hinges-locked 
vehicle modal data for the spacecraft were 
obtained by using NASTRAN. From a detailed finite 
element model having 19,434 degrees of freedom and 
3,239 nodes, 63 hinges-free and 67 hinges-locked 
elastic modes below 25 Hz were computed. Since the 
vehicle is essentially symmetric (the sensor 
causes a slight asymmetry), both symmetric and 
ant i symmet ric vehicle modes arise in transverse 
bending and in-plane bending of the arrays, and 
the vehicle modes are categorized accordingly 
in Table 1 and Table 2. Fig. 3a confirms the 
prediction from the identity (ill) that the 
hinges-free modal coefficients, in this case 
♦ ^ (u=l , . • • ,63) for the y^-solar array, form a 
nonconverging series. In Fig. 3a, the largest 
modal coefficients for p-8 , 11 , 18 ,28, . . . cor- 

respond to those vehicle modes which predominately 
entail torsion of the y^-array about y^-axis 
(Table 1). In contrast, those contributing to 
0^, namely, , . . . , 63 ) , form essentially a 

convergent series because the sensor is rigid, and 
symmetric transverse bending of the arrays 
(Table 1) or local high-frequency deformation of 
Bq at the sensor base produce $ (p=l , . . . ,63) . 
The hinges-locked coupling coefficients h, 
for y^-array and h^_ for 0 


rotation for^tKe 
Eig • 4. 
converging series. 


3ay 

modes a 3 l,...,67 are 

Unlike ♦. , forms a 

lpy lay 


. 3 X 

displayed in 


The identities (III) and (V) are the simplest, 
for they involve only modaL coefficients, no 
frequencies. The identity (III) is illustrated in 

HF 

Fig. 5. The error indexes e^ k (k-I, 3) (HF means 

hinges-free) are the corresponding diagonal 

elements of the (4x4) matrix [l *3 $T<t> l” 1 . In 

— — —A— A 

contrast to their zero ideal vaLue, the asymptotes 


EARTH 

o 



Figure 2. A Four-Body Deformable Spacecraft 
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Table 1. Hinges-Free Modes 


Characteristics 
of the Mode 

Mode No, 

Affected 
Rotational 
Degrees of 
Freedom 

Transverse bending of 
arrays and A-frames 

Symmetric 

1. 5,9, 13, 15,21,... 

0oy' 03 y 

Antisymmetric 

2, 6, 10, 14, 16, 22, ... 

0o z 

Torsion of arrays and 
A-frames 

Array 1 

8, 11, 18,28,... 

0 1 y 

Array 2 

7, 12, 19, 29,... 

02y 

In-plane bending of 
the A-frames and 
solar arrays 

Symmetric 

3,... 

None 

Antisymmetric 

4, 17, 20,... 

0ox. ©3x 


Table 2. Hinges-Locked Modes 


Characteristics 
of the Mode 

Mode No. 

Affected 
Rotational 
Degrees of 
Freedom t 

Transverse bending of 
arrays and A-Frames 

Symmetric 

1,2, 6,8, 12, 16 

0oy 


1,2,... 

©3y 


Antisymmetric 

3, 9, 13, 17, 20, 27,... 

0OZ 

Torsion 

Array 1 

5,6, 10, 11, 14, 15,22,23,... 

e 1 y 

In-plane bending of 
arrays and A-frames 

Antisymmetric 

7, 18,21,... 



18,... 

©3y 

Vibrations of the spacecraft 

28, 35, 36, 37, 41 , 42, 43, 44, 
47, 48, 49, 53, ... 

©ox 



28, 36, 37, 41 , 42, 43, 44, 45, 
47,48,49, 53,... 

0oy 



28,35, 36.37.38,41,42,43, 
44,47,48, 53,... 

0oz 

t Information about the interaction with 02y and 03 x not available 


W 





i+3„yl 



Figure 3. Hinges-Free Modal Coefficients of Articulation Motion of y j- Solar Array, 
blyyl» and of Sen » or B 3* l*3yyl > about y 3 -Axi* 
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or * MODE NO. 

Figure 4. Hinges-Locked Modal Coefficients 
Associated with the Articulation Motion 



fi * MODE NO. 


Figure 5. Hinges-Free Identify III: Diminishing 
of Error Index with Hinges-Free Modes m 

of these indexes are 0.0406 for k»l (8. rotation 

iy 

of the y^-array) and 0.689 for k=3 (8-^, the 

x-rotation of the sensor). The error indexes 

diminish discretely at appropriate modes as 

predicted by Table 1. For instance, for 8, , the 
HF . ly 

error index e^ diminishes at the torsional modes 

p a 8, 11 , 18,28, ... . The index for the y 2 “array 

motion (k=2) is the same as that for k-1 , except 
that it decreases instead at the adjacent 

torsional modes p=7 , 12 , 19 , 29 , . . . (see Table 1). 
Surprisingly, the asymptote of the error index 
for 8^ (k=4) , not included in Fig. 5, hovers at 

0.9976 instead of decreasing to the ideal value 

zero. Fig. 6 illustrates the hinges-locked 

T 

cl cl -1 

identity (V), rearranged as h A ^ * _1. For 

discussing this identity and the ones follow- 
ing, define a "completeness index 0 11 which 
approaches unity for an error-free model [Refer- 
ence 1]. [The completeness index for Fig. 5 is 



0 1 10 20 30 40 


a -MODE NO. 

Figure 6. Hinges-Locked Identify V: Crowth of 
Completeness Index with Hinges-Locked Modes a 


(1 - e^)*! The growth of the diagona^l elements 

(1,1) and (4,4) of the matrix *are 

depicted in Fig. 6, Surprisingly, alongside the 
HF HL 

error index in Fig. 5a,C^(HL means hinges- 

locked) approaches unity in just two hinges-locked 
torsional modes, 3 and 6 (Table 2), and its 

asymptotic value is 1.029. Furthermore, by con- 
trast with the hinges-free completeness index 
HF 

C/# equal to 0.0024 (that is, the above mentioned 
44 HF 

error index e / of 0.9976), the hi nges-locked 
44 

completeness index 35 0.9421 in Fig. 6b is 

remarkable; in fact, the first hinges-locked mode, 
a symmetric transverse bending mode of the arrays 
(Table 2), contributes a mighty share, 0.9412, 


The identities which involve frequencies as 
well are now illustrated. First, consider the 
identity (VII) which is summed over ail hinges- 

free modes (p=l , . . . , 63 ) for a specific w . When 
c u 

u and are the same to several decimal 

peaces, it is difficult to verify this identity in 
this form. On the other hand, the identity (VIII) 
indicates that when w and w C are truly the 

U a 

same, the corresponding poles and zeros cancel 
each other without affecting the articulation 
dynamics. A physical explanation of this is that 
when hinges-free and hinges-locked frequencies are 
truly equal, that particular mode does not 
contribute to the articulation motion, so such a 
mode may be deleted from the study. In numerical 
work, however, it is difficult to establish true 
equality between two real numbers. Besides, as 
will be seen shortly, for the example in hand, 

sometimes even though w and w C are the same up 

V ci 

to three or four decimal places, the minuscule 
difference between the two is still important for 
the verification of an identity. Consequently, the 
following results are obtained without truncating 
either modal set. Returning to the identity VII, 
one finds that when a > 9, hinges-locked 

frequencies w C are so close to a corresponding 
hinges-free frequency u that the determinant, 
instead of being zero, becomes an arbitrarily 
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large number. Among a=l,..., 9, Che identity (VII) 
is best satisfied with a=7 and next best with 
a=2 , for wnich the determinants are, respectively, 
-0.00415 and 0.04732 (Table 3). The circumstances 
which produce these results are revealed by the 
identity (X). For a given u C , when all available 
hinges-free modes are added to calculate the 
U,k) element of the left side of (X), it is 
HF 

denoted C where "asy means asymptotic 

4k, asy HF 

value. Fig. 7 shows a for fc*l*3, and 4, The 

ideal value of this index is unity; however, 
when w =* (j C for some u and a, this index 

assumes an arbitrarily large vaLue, and for 
plotting purposes, such large numbers are replaced 
by 2 without aLtering their signs. In the 
left, side of Fig. 7. a, in the useful range 0 


to 1, the maximum value of concerning the 

yj-array rotation, 0^ , is 0.451 for the hinges* 
locked mode a-2 — a symmetric transverse bending 
mode of the arrays (Table 2). On the other hand, 
the first torsional hinges-locked mode having 
significant coupling with the rotation 6. is 

HF *y 

a=5 (Table 2 and Fig. 4), but C ^ a correspond- 
ing to a— 5 is 0.16, less than 0.451 for a=2. 
HF 

Although the index C ^ agy for a “ 5 should be, 

intuitively, greater than that for a=2, this 
does not happen because the hinges-free frequency 
w L (0.25724 Hz) is close to (0.25719 Hz). To 

determine the contribution of the mode y = l, the 
growth of Cjj with successive addition of y to the 
asymptotic value 0,451 for a=2 is shown in the 
right side of Fig. 7. a. C ^ is found to escalate 


discretely at y = l ,8 , 11 , 18, 28 ,29,35 ,41 ,42 , . . . , 
which, except for y=l, involve torsion of the 
array 1 (Table 1). The contribution from the 
hinges-free mode y=l, a symmetric transverse 
bending mode of the arrays (Table 1) like 
a=2 hinges-locked mode, is however, extraordi- 
narily large: 93X. Nevertheless, the bending 
mode u*l is not pertinent to the articulation 


motion 8. » soQn a 0.451 for a*2 cannot be 

ly' 11 

accepted, and, instead, ■ 0.16 for o*5, a 

torsional mode is accepted. Next consider the 
sensor motion 0_ . The corresponding index, 
HF 

(L, , shown on the left side of Fig. 7.b, is 

33, asy' HF 

1.0059 for 0*7 (compare with C ^ agy > and recall 

from Table 3 the value 0.00415 of the identity VII 

HF 

for a=7). The growth of versus y for a=7 is 

displayed on the right side of Fig, 7.b, where it 

is observed to become unity at once when y=4. To 

understand this, note that both y=4 and o=7 modes 
involve antisymmetric in-plane bending of the 
arrays — a motion which induces 0^ (see 
Table 1 and Table 2), and that w, * 0.59541 and 
* 0.59538 Hz. When the rotation 9^ x of the 

rigid sensor is locked, the moment of inertia 
which must be turned by the antisymmetric in-plane 
bending, is increased, and that lowers the fre- 
quency commensurately . The ratio of the moment of 
inertia of the sensor and of the core body, both 
about XQ-axis, is 0.0717. The decrement of 3.0E-5 
Hz noted above in the frequency u>, is mathemat- 

HF ^ 

ically so precise that C ^ becomefl unity at once 
when a=7. Moreover, although and are the 


same up to tnree decimal places, the two modes 
cannot be truncated from the study of the 
verification of the identities VII and X, Next, 
consider 8^ rotation of the sensor — the rotation 
coupled wi t¥i the transverse symmetric bending of 
the arrays (Table 1 and Table 2). The associated 

Hp 

index. , versus a is shown in Fig. 7c. In 

9 v '44,asy’ 

the range 0 to 1, the most it becomes is a 
startling low value: 0.07836 for a*2; for this a, 
HF 

the growth ofC^ with y indicates that 99.99X 

contribution arises from the first symmetric 
transverse bending mode y=l. 


The verification of the identity (IX) for 
l=k=l and 4 is considered in Fig. 3. Since this 
identity relates to hinges-locked modaL para- 

hl 

meters, its left side is denoted Earlier, the 

identity (V) and Fig. 4 established that the 
hinges-locked coupling coefficients form a con- 
verging series. Therefore, the determinant 
HL 

identity (VI) and C ^ asy Fig. 8 do not 

become arbitrarily large numbers once y > 28. 

Indeed, only for y=3,4,21 ,26,27, is the index 

unbounded, by contrast with the hinges- 
^ik,asy up 

free index Ci i in Fig. 7a which is unbounded 
11, asy Hl 

for all a >9. The index (J depends on the 


selected hinges-free 
having 


frequency u> » for a' 9 
2W 2, 


u» C >w , the term f 1 - w C /w") becomes 
ay' 'ay 

negative and these particular hinges-locked 

modes diminish the sum. Focusing first on 

HL 

C, , , surprisingly, it stabilizes early on to 

w ll,asy r 

1.05 when y=7 or 8 — the first two hinges-free 

H L 

torsional modes. The ascent of to 1.05 for 

y*8 with hinges-locked modes a (Fig. 8a) indicates 
significant contributions from a*5,6,10, and 
11 — all torsional modes (Table 2); the con- 
tribution from higher torsional modes 

attenuates rapidly because of the fast convergence 

of h? 1 . As for the rotation 8, , the maximum 

-lay lit . v 

value of C // i displayed in Fig. 8b, in the 
44, asy 

range 0 to 1 is 0.959 when u=5 — the second hinges- 
free symmetric transverse bending mode of the 

HL 

arrays (Table 1). The growth pattern of versus 

a for y = 5, also shown in Fig, 8b, states chat 
virtually the entire contribution arises from the 
first hinges-locked mode (a*l) involving symmetric 
transverse bending of the arrays, 

VI. SUMMINC UP 

To draw conclusions about the relative merits 
of hinges-free and hinges-locked vehicle modes, 
Table 4 summarizes the completeness indexes for 
the identities (III), (V), (IX), and (X), 

Evidently, the hinges-locked indexes are far 
closer to unity than the hinges-free indexes. The 
superiority of the hinges-locked vehicle modes to 
the hinges-free modes is established most 

persuasively by comparing the indexes for the 
articulation motion 6_ of the sensor: 

HF 3y 

Cm. are 0.0024 and 0.0784 — far remote from 
44, asy ut 

unity, whereas aay ere 0.9421 and 0.9593 — 

almost unity. It must be understood, nevertheless, 
that the identities (X) and (IX) (or VII and VI) 
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table 3. Identity VII: Variation of the Hinges-Free Determinant With Hinges-Locked 

Modes; IdeaL Value » U 

« 1 2 3 4 5 6 7 8 9 

VII 0.97105 0.04732 0.27012 0.79457 0.68305 0.68009 -0.00415 0.19551 -0.32426 

Minimum value of the determinant among those for a = 1 0, .... 63, is 74.4, and the maximum value is 
oo when Wp = «£ up to several decimal places 



HL MODE NO. a HINGES-FREE MODES 



10 20 0 10 20 30 40 i 

HL MODE NO. a HINQES-FREE MODES 


Figure 7. Identity X: Asymptotic Values of the Hinges-Free (HF) 

Completeness Indexes Versus Hinges-Locked (HL) Mode a, 
and Growth of this Index Versus Hinges-Free Modes 



HINGES-FREE MODE NO.^ 



HINGES-FREE MODE NO. n 



RETAINED HINGES LOCKED MODES 


Figure 8* Identity IX: Asymptotic Values of the Hinges-Locked 

Completeness Indexes Versus Hinges-Free Modes, and 
Growth of this Index with Hinges-Locked Modes 
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Table 4. A Summary of the Completeness Indexes for 
Hinges-Free and Hlnges-Locked Vehicle Modes; 
Ideal Value - 1 


Hinges- 
Free (HF) 
Indexes 

Identity 

III 

1 _ e HF 
1 0 U . asy 

Identity 

X 

Hinges- 
Locked (HL) 
Indexes 

Identity 

V 

Identity 

IX 

Articulation 

Motion 

Associated 
Mode of 
Deformation 

° 1 1 ,asy 

0.9594 

0.160 

rH L 

u 11, asy 

1.029 

1.05 

© 1 y 

Torsion 

(HF 

u 33,asy 

0.311 

1.0059 

^*"33, asy 

TBDt 

TBDt 

®3x 

Antisymmetric 
in-plane bending 

(HF 

44,asy 

0.0024 

0.0784 

(f L ,„ 

u 44, asy 

0.9421 

0.9593 

e 3v 

Symmetric trans- 
verse bending 

|to be deterr 

nined 


represent two different situations: in the former, 
the hinges-free modes are employed to yieLd a 
bounded response at a hinges-locked frequency; and 
in the latter, the hinges-locked modes are used to 
elicit an unbounded response at a hinges-free 
frequency. Therefore, a comparison of the indexes 
from these identities is slightly inappropriate 
perhaps; yet the conclusion from Table 2 seems 
inevitable that the hinges-locked vehicle modes 
yield a much more accurate model for simulation 
than the hinges-free vehicle modes do. This is 
caused by the nonconvergence of the hinges-free 
modal coefficients in contrast with the rapid 
convergence of the hinges-locked coupling 
coef f icients — the attributes corroborated by the 
identities. Besides contrasting one family of 
modes with the other, the identities are clearly 
useful in sifting through scores of finite-element 
generated modes to select a few pertinent modes 
for an articulation degree of freedom in consid- 
eration. An important extension of the preceding 
work is to devise identities which involve modal 
coefficients and frequencies of only one family of 
modes, hinges-free or hinges-locked, not both. 
Hughes^ has formulated such identities for an 
elastic body with no articulated members. 
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